A simple model for the vibrational modes in honeycomb lattices 
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The classical lattice dynamics of honeycomb lattices is studied in the harmonic approximation. 
Interactions between nearest neighbors are represented by springs connecting them. A short and 
necessary introduction of the lattice structure is presented. The dynamical matrix of the vibrational 
modes is then derived, and its eigenvalue problem is solved analytically. The solution may provide 
deeper insight into the nature of the vibrational modes. Numerical results for the vibrational 
frequencies are presented. To show that how effective our method used for the case of honeycomb 
lattice is, we also apply it to triangular and square lattice structures. A few suggested problems are 
hsted in the concluding section. 
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I. INTRODUCTION 

In crystals, atoms vibrate about their equilibrium po- 
sition. Calculation of the vibrational frequencies and 
modes is an important and presently actively studied sub- 
ject in solid-state physics. To interpret various properties 
of crystal lattices — for example, specific heat, thermal 
expansion coefficients, elastic constants — it is essential to 
take into account the lattice vibrations. Fortunately, sev- 
eral good textbooks of the subject are available (here we 
give just a small selection of the vast literatur o^i^i'^i'^i^ ), 
in which the quantum and classical treatments of lat- 
tice dynamics are presented lucidly. The vibration of the 
atoms depends on the interatomic interaction within the 
crystal. The interatomic interaction potential is often 
approximated by including only quadratic terms of the 
displacement of the atoms. This is called the harmonic 
approximation, and is the usual starting point for devel- 
oping the theory of lattice dynamics. In the framework 
of classical mechanics, lattice dynamics can be treated 
through the equations of motion for the atoms. To de- 
termine the vibrational frequencies and the correspond- 
ing modes one needs to calculate the eigenvalues and the 
eigenvectors of the so-called dynamical matrix, which can 
be obtained from the interatomic interaction potential. 
If the dynamical matrix is known, the eigenvalue prob- 
lem is straightforward, though only numerical solutions 
are available in many cases. Numerous textsiiSiSi^iSiSil 
provide problems for the calculation of the vibrational 
modes in different crystal lattices. Fortunately, there 
are some well-known examples — such as linear chains 
with/without bases, square and cubic lattices — in which 
analytical methods can be usedi*2iiii^. These exam- 
ples are important for a deeper understanding of the the- 
ory of lattice dynamics. 

In this paper we shall give another nontrivial example, 
namely the lattice dynamics of the honeycomb structure 
shown in Fig. ^ The classical lattice dynamics of the 
honeycomb lattice is modelled by using the harmonic ap- 
proximation and assuming only nearest neighbor inter- 
actions between the atoms. In this simplest model the 
interaction between nearest neighbors can be represented 



by springs connecting them. It is assumed that black and 
white atoms have different masses, mi and m2, respec- 
tively. The problem can be solved analytically and its 
solutions may provide insight into the nature of the vi- 
brational modes. As it will be seen below, the analytical 
method is possible because of the symmetry properties of 
the honeycomb lattice. As a pre-requisite to the solution, 
a short and necessary introduction of the lattice struc- 
ture is presented. More details on crystal structures can 
be found in many books on solid-state physicsiiSiSili^. 
Finally, in connection with the honeycomb lattice, the 
intensive research on nanotubes^ should be mentioned. 
Nanotubes are quasi-one-dimensional cylindrical struc- 
tures in which a two-dimensional honeycomb lattice of 
crystalline graphite is rolled up into a cylinder. The 
work presented in this paper may serve as an introduc- 
tory study of the far more complex lattice dynamics of 
graphite and nanotubes. 




FIG. 1: Honeycomb lattice with two different atoms (black 
and white circles). The straight lines correspond to un- 
strained springs of force constant D connecting the nearest 
neighbor atoms in equilibrium state. The masses of the black 
and white atoms are mi and m2. 

The general method developed for calculating the vi- 
bration modes of honeycomb lattices proves to be use- 
ful for studying other lattice structures. Therefore, we 
shall present a few numerical results for the vibrational 
frequencies of equilateral triangular lattices and square 
lattices in which both the nearest neighbor and second 
nearest neighbor interactions are taken into account. 

Our approach used in this paper may provide a didac- 
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tically useful examples for introducing the basic lattice 
dynamical notions developed in solid-state physics. 

The rest of the text is organized as follows. In Sec. 
nil a formula for the spring force acting on the atoms is 
derived. In Scc. lIIII the description of the honeycomb lat- 
tice and the equations of motion for the atoms are given. 
In Sec. lIVI our analytic method for solving the eigenvalue 
problem is presented. In Sec. IVl some numerical results 
are mentioned. A few results are presented in Sec. |^ 
for the case of triangular and square lattices. The con- 
clusions and some suggested problems are given in Sec. 

Ivlil 



vector u, and that the constant of proportionality is the 
product of the spring constant D and the direct product 
of the unit vector d/|d| with itself. In the following sec- 
tion the equations of motion of the atoms are formulated 
by using the above form of the spring force acting on the 
atoms. 



III. THE EQUATIONS OF MOTION 

In Fig. n the position vector R of the black atoms in 
the honeycomb structure can be given in the form 



II. THE FORCE 

First, the force acting on an atom by a spring connect- 
ing the atoms is calculated. The springs are assumed to 
be ideal, ie. the spring force is proportional to the change 
in the spring's length. Denoting the vector between the 
two endpoints of the unstrained (strained) spring by d 
(r), the change in the length of the spring is |r| — |d|. 
If the spring constant of the ideal springs is D then the 
restoring force is 



F = -D(|r|-|d|) 



(1) 



It is useful to express the force in terms of the dis- 
placement vector u = r — d. In the harmonic ap- 
proximation the strain of the spring is assumed to be 
small, ie. second order terms are neglected in u. Ex- 
plicitely: |u|/d — |r — d|/d ^ 1. Thus, we have |r| — 



du^ 
"32- J 



and the force can be writ- 



ten as 



-D 



(du)d 



(2) 



It is convenient to define the outer or direct product a®b 
of two vectors a and bifl*ii. This is a matrix whose a, (3 
element is a product of the ath component of a and the 
/3th component of b, that is 



(a(g)b; 



Q/3 



a* 6/3, 



(3) 



where * denotes the complex conjugation (of complex 
vectors). The following direct-product identities will 
prove useful in the calculations below: 

(a(g)b)c = a (be), (4a) 
(a(g) b)(c (g) d) = (be) a (Kid. (4b) 

In direct-product notation, the ^ force can be rewritten 
as 



(5) 



From this form it is clearly seen that the force in the har- 
monic approximation is proportional to the displacement 



R 



(6) 



where ai and a2 are the independent primitive trans- 
lation vectors shown in Fig. [21 while ji and j2 range 
through all integer values (i.e. positive and negative in- 
tegers, as well as zero). The primitive translation vec- 
tors have the same magnitude |ai| = |a2| — \/3a, where 
a — |di| is the length of hexagons' side. The primitive 
cell of the honeycomb lattice can be chosen as the paral- 
lelogram whose non-parallel sides are the two primitive 
translation vectors ai and a2 (in Fig. |21the sides of the 
parallelogram are ai , a2 and the two dotted lines) . In the 
honeycomb lattice each primitive cell with a two-point ba- 
sis contains two atoms (black and white circles in Fig.[T]l: 
one is at the common starting point of vectors ai and a2 , 
while the other one is at 1/3 (ai + a2) in Fig. |21 




FIG. 2: The primitive translation vectors of the honeycomb 
lattice are vectors ai and a2; the primitive cell is a parallel- 
ogram whose sides are these two vectors and the two dotted 
lines. The three vectors di,d2,d3 will be useful in the calcu- 
lations below. 



by 



It is useful to introduce three unit vectors (rii, n2, na) 



n, = 



i = 1,2,3. 



(7) 



The vectors di, d2, da, are shown in Fig. [2 their magni- 
tude is |di| = a. Thus, the primitive translation vectors 
are ai = d2 — di = a(n2 — rii) and a2 = ds — di = 
a(n3 — rii). Later on the following identities for vector n 
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will frequently be used: and E is the 2x2 unit matrix. 

3 

,g , The displacements of the two atoms in the primitive 

~ ' ^ ' cell specified by the lattice vector R are denoted by u(R) 

andv(R). From lO it can be shown that the components 
3 of the displacements perpendicular to the plane do not 

n 6?) n = — F (f^h) J' J' 

iij iij ^ ji;, \^ou) contribute to the force in the harmonic approximation. 

Each atom interacts with its three neighbors. Denoting 

{1 if J = k, the masses of the two basis atoms by mi and m2, and 

^ , (8c) using ISJ for the spring force acting on the atoms, the 

^2 j ^ following equations of motion are obtained for u(R) and 

v(R): 



E 



mi u(R) = -£> ni ni [u(R) - v(R)] - Z? na ® na [u(R) - v(R - ai)] - L> ng ® ng [u(R) - v(R - as)] ,(9a) 
m2V(R) = -Dni^ni [v(R) - u(R)] -£in2(8)n2 [v(R) - u(R + ai)] -Dna^ng [v(R) - u(R + as)] .(9b) 



Eq. Q is an infinite set of equations for the displacements 
u(R) and v(R). Following the traditional method, we 
seek a solution representing a wave of angular frequency 
Lu (q) and wave vector q: 



where 



u(R) = liM e-(q)*+^qi^, 

V(R) = ^ (q) ^ii^(q)t+iqR^ 



(10a) 
(10b) 



where the vectors u (q) and v (q) are to be deter- 
mined, and q is in the first Brillouin zone. The 
usual Born-von Karman periodic boundary conditions 
are applied^^i?i'^i'?i?, i.e. u(R + iV,aj) = u(R) and v(R + 
NiSLi) = v(R), where A^i and N2 are large integers whose 
product N1N2 ^ N is the total number of primitive cells 
within the crystal. The periodic boundary conditions re- 
strict the allowed wave vectors q in Eq. 1)10(1 to the form: 



^^iVi^^ + iV2^^' 



(11) 



where pi and p2 are integers, and the are the reciprocal 
lattice vectors defined by a^bj = 27r5y, i,j = 1,2. It is 
convenient to choose pi , P2 such that q is limited to the 
first Brillouin zone. For example, if Ni and N2 are even 
integers (which will be irrelevant in the limit N — > 00), 
-NJ2 < Pt < Ni/2 {i = 1,2). More details on the 
concept of the Brillouin zone can be found in many books 
on solid-state physicsiiSiSi^i^. 

Substitution of the solutions (|10|l into the equations of 
motion ^ leads to an eigenvalue problem: 



D(q) 



u(q)' 
v(q) 



= ^2 (q) 



u(q)' 
v(q) 



(12) 



D(q) = 



Dii(q) Di2(q)' 
D2i(q) D22(q) 

D ^ 

Dii(q) = y^nj®nj = 

nil ^-^ 

D ^ 

D22(q) - — E"J®"J' = 



mi — ' ZTOi 



m2 



Di2(q) 
D2i(q) 



D 



3D 
2m2 



-E, 



(13a) 

(13b) 
(13c) 



Y^mlm^ 
[Di2(q)]*. 



(13e) 



In the last equality of (|13b|l and H13c|l we used Eq. I|8bp . 
The matrix D(q) is commonly called the dynamical ma- 
trix. The eigenvectors of the dynamical matrix are the 
unknown vectors u (q) and v (q) of Eq. (|10|) . In our case, 
D(q) is a 4 X 4 matrix, while the (q) are 2x2 matrices 
(i,j = 1,2). For a given q there are 4 eigenfrequencies 
and 4 corresponding, mutually orthogonal eigenvectors, 
which give the amplitudes of the wave solutions l|10|l . 
In the literature these eigenvectors are called vibrational 
modes. 

The eigenfrequencies can be found by solving the fol- 
lowing equation for A(q) ~ uj^ (q): 



det [D(q)-A(q)I]=0, 



(14) 



where I is the 4x4 unit matrix. The 4x4 determinant 
can be evaluated directly; this is, however, rather tedious. 
Below we present another approach, which allows us to 
find analytical solutions for the eigenfrequencies. Note 
that the dynamical matrix is hermitian (see Eq. (|13e|l l. 



4 



therefore its eigenvalues are real. The basic idea is that 
D(q) is a block matrix, whose elements are the 2x2 
matrices Dy (q), and Dii(q) and D22(q) are diagonal. 
This way, the original eigenvalue problem is reduced to 
finding the eigenvalues of a 2 x 2 matrix. This problem 
can be solved analytically. 



IV. ANALYTIC SOLUTION OF THE 
EIGENVALUE PROBLEM 



where 

ri{q) = cosqai + cosqa2 + cosq(ai - a2), (20) 

has been expressed in terms of the primitive translation 
vectors ai and a2. Thus, the eigenvalue Ai corresponding 
to the eigenvector vi is 



Ai 



3 + 27;(q) 



(21) 



In this section an analytic method is presented to de- 
termine the eigenfrequencies and the vibrational modes 
of the honeycomb structure. Since Dii(q) and D22(q) 
are diagonal matrices, one can eliminate u(q) from 
Eq. lO, and find by using Eqs. lfT3B)l - lfnc|l 



PP*v(q) = A(q)v(q), where 



(15a) 
(15b) 



Eq. H15a|) is another eigenvalue equation. The matrix 
PP* can be rewritten as 

3 

PP* = [(n,®n,)(nfe®nfc)]e''i(<*^-'^'=) 

j,k=l 
3 

= ^ [(n,®nfe)(n,nfe)]e^'i('i^-'i'=) 

3 ^ 3 

3 = 1 j,k=l 
33 ^3 

i=i j,k=i 



9 1 

-E - - f ® f, 

4 2 ' 



where 



(16) 



(17) 



In the derivation, identities (j4b|l . (|8b|) and (j5c|l have been 
used. Note that f is a complex vector. 

From H4a|l it is easy to see that vi = f is an eigenvector 
of the matrix PP*: 

PP*Vi = PP*f = ^ f - i f2 f = Q - i VI, (18) 

with eigenvalue Ai — 9/4— 1/2 f^. Using the dot 
product can easily be evaluated: 

3 , 3 



3,k=l 

= 3-77(q), 



j.k=l 



(19) 



The other eigenvector must be perpendicular to vi. 
It is useful to introduce three unit vectors, Ii,l2,l3 (see 
Fig. El, such that 



XjXij = 0, 
Ij-nfcl = if j ^ k. 



(22a) 
(22b) 

(22c) 



Next, the vector 




FIG. 3: The three vectors Ii,l2,l3 are obtained by a 90° 
rotation from ni, n2, na. 



(23) 



is defined. Note the presence of the vectors dj in the 
exponent. Then, from (|22|l it is obvious that g is per- 
pendicular to the vector f: 

3 

f g = ^ njhe-''^'^^-''^'^'' = 0. (24) 

3,k=l 

Therefore, the other eigenvector of PP* can be chosen 
as V2 = g, and the corresponding eigenvalue A2 can be 
found from PP*V2 = PP*g = |g- 5f(fg) = j^2- 
Thus, the second eigenvalue is 



A2 



9 



(25) 



Note that this eigenvalue is independent of q. Finally, 
using (|15c|l one can obtain the four eigenfrequencies of 
the dynamical matrix D(q): 



5 



Ai,2(q) 



1 1 
1 

mi m2 



± 



1 1 

mi m2 



12 + 8 77(q) 

9 17111712 



(26a) 



W3(q) = V'^3(q) = and W4(q) = \/Al(q) = W-Z? h , 

y 2 V mi m2 / 



1 1 



(26b) 



where Ai and A2 (± signs, respectively) correspond to Ai, 
while A3 and A4 to A2. 77 (q) is defined in if^ . 

To find u(q) one can use Eqs. ifT^ and Ijl3b|l - H13c|l . 
and obtain 



"(q) 



D/ ^771 im2 

'A(q) - 3i:>/(2mi) 



'q<iip*v(q). (27) 



First, P*vi and P*V2 are calculated. Using |ga|l, (jHa] 
and l(5c|l 



P*vi = P*f = ^nj(8)nje-'i<^^^nfc( 



i=i 



fc=i 

„jq(dfc-d3) 



i=l j,fc=l 
3 3 



fe=l 



where c = — ^ X]fe=i 6*'''*'° is a scalar value depending on 
q. Similarly, 



P*V2 = P*g = ^nj®n,e~''J'*^^l 



fee 



fe=i 



--^ (l.e^qdi +l2e»<id. ^I^g.qd3)^_3^* (29) 

where we have made use of the identities n2 — rii = a/S I3, 
n3 - n2 = VSli, ni - n3 = ^312 and di + d2 + d3 = 0. 
Finally, using Eqs. (|27|l - 129|) the four eigenvectors are 



UJ3 



-cCi(q) P 
f 



'"1 „*giqdi 



W2 



-cC2(q) r 
f 



/ "12 „* iqdi 

Y 7771 S ^ 



Where Q (q) ^ , ^/V™ ^-'^''^ 
'^^^ A, -3i?/(2mi) 



(30a) 



(30b) 



(30c) 



and f , g, Ai, and c are functions of q defined in Eqs. (|17|l . 
(1^ . if^ and after Eq. (j^Hl, respectively. It can be 



shown that the four eigenvectors are mutually orthogo- 
nal. 

It is also worth mentioning that the eigenmode cor- 
responding to ijJ3(q) = is a special one, for example, 
[u(R) — v(R)] rii = for this mode. This can easily 
be proved by using pO|) and the expression for the UJ3 
eigenmode given by (|30b|) . This means that neglecting 
the second order terms in the displacements u(R) and 
v(R) all springs are unstrained, ie. their lengths remain 
a. The atoms of the whole lattice move as if the spring 
system were a universal joint. This movement of the 
atoms does not require any energy (the eigenfrequency 
is zero), therefore there is no resistance to a shear of 
the lattice. Consequently, the stability of the lattice is 
lost, and the whole honeycomb structure becomes unsta- 
ble. In nature, besides the central interaction between 
nearest neighbors there exist tt,*'* neighbor (77, = 2, 3, • • • ) 
and lateral interactions as well'*'^. In the latter case, the 
interatomic interaction energy depends on the angle be- 
tween two adjacent bonds. Such interactions give rise to 
the so-called bond-bending forces. These additional in- 
teractions stabilize the lattice and one must take them 
into account to calculate the vibrational modes more ac- 
curately. Nanotubes provide a good example for this^. 



V. NUMERICAL RESULTS 

In this section a few numerical results for the eigen- 
frequencies are presented. As it is seen from (|26|) only 
ti;i(q) andw2(q) have nontrivial q-dependence. In Figs. 01 
and O the contour plot of wi(q) and W2(q) are shown for 
mi/m2 = 2. The vector q = {qx, Qy) is taken in a Carte- 
sian coordinate system in which the x-axis is parallel to 
the direction of the vector — di. In these contour plots the 
lines show the constant eigenfrequencies in the parameter 
space q = [qx,qy)- It is seen from the figures that a;i(q) 
has a deep minimum (zero) at q = 0, while 1^2 (q) bas 
a maximum there, equal to 104^ = y^3D{l/mi -\- l/m2)/2 
for arbitrary mass ratio mi/m2. The contour plot is also 
a useful representation of the vibrational modes to show 
the symmetries of the eigenfrequencies in the q space. 
One can see from Figs. 0] and |5l that both wi(q) and 
0-12 (q) have, eg. a 60° rotation around the axis which 
goes through the point q = and perpendicular to the 
{qx,qv) plane. 
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FIG. 4: The contour plot of lui (q) (in units of ^ D/m\) versus 
q (in units of 1/a) for m\/m2 — 2. Qx is parallel to the vector 
-di. 



In Fig. El the eigenfrequencies ivi (q) , uj2 (q) and (q) 
are plotted along the lines in the first Brillouin zone join- 
ing the points T, M and K shown in the inset of the 
figure. In solid-state physics it is a common practice to 
plot the eigenfrequencies along the lines between high 
symmetric points. Owing to the symmetry properties of 
the eigenfrequency in the q space it is enough to calcu- 
late the eigenfrequencies only inside the area enlosed by 
these lines. In the same Cartesian coordinate system as 
in Fig.mthe points T, M and K are qr = 0, qw = 0) 
and qK = (|f , g^^ ), respectively. Here a is the length 
of hexagons' side. Note that W4(q) is independent of q. 

In the literature the vibrational mode corresponding to 
u;i(q) is called the acoustical branch because the eigen- 
frequency at small |q| has the same form as that of the 
sound waves, namely uj = u|q|, where v is the sound ve- 
locity. Those vibrational modes for which the frequencies 
does not tend to zero as |q| — > form the optical branch. 
In our example a;2(q) and Ci^4(q) belong to the optical 
branch. It is clear from Fig. that at the point K there 
is a frequency gap between the acoustical and optical 

branch. This gap is equal to ^ 1 1 / y/nii — 1/ | . 

We now take mi = m2 corresponding to a single 
graphite layer. A similar plot as in Fig. is shown in 



FIG. 5: The same contour plot as in Fig. |l|for aj2(q). 
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FIG. 6: The eigenfrequencies (in units of \J D lm\) lo\(c\) 
(solid), a'2(q) (dashed), W4(q) (dot-dashed) along lines be- 
tween points r, M and K. These points are in the first Bril- 
louin zone shown in the inset. The mass ratio is mi/m-z — 2. 



Fig.[7|for this case. One can see that the gap disappears. 
These results are similar to those obtained from the re- 
cent accurate calculations for two modes of the isolated 
graphite layer—. 
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FIG. 7: The same plot as in Fig. |S|for m\/m2 — 1. 

VI. OTHER LATTICE STRUCTURES 

To demonstrate the effectivity of the general method 
developed in Sections ITll and [llll we now present further 
examples of lattice structures. 

In Fig.|Slan equilateral triangular lattice is shown. We 
assume only nearest neighbor interactions repersented 
by springs of force constant D. One can easily write 




FIG. 8: Equilateral triangular lattice with nearest neighbor 
interactions. The straight solid lines correspond to unstrained 
springs of force constant D connecting nearest neighbor atoms 
in equilibrium state. The unit vectors ni , 112 and ns corre- 
spond to the directions of the unstrained springs of length a. 



down the equations of motion for the displacements of 
the atoms using Eq. ||SJ) and it yields the 2x2 dynamical 
matrix: 

2D . 

D(q) = [1 - cos(qnia)] rii (g) rij, (31) 

m ^ — ' 

i=l 

where the unit vectors (z = 1, 2, 3) are shown in Fig.|Sl 
The eigenfrequencies are again the square-root of the 
eigenvalues of the dynamical matrix as in (jl4(l . From a 
simple calculation one finds 



\ 



Dn + D22 ± {Dii - D; 



^22) 



(32) 

where D^j is the ij matrix element of the 2x2 dynamical 
matrix D(q) for a given wave vector q. In Fig. |51 the 
eigenfrequencies wi(q) and W2(q) are plotted along the 



lines in the first Brillouin zone joining the points F, K 
and M shown in the figure. In the Cartesian coordinate 
system in which Qx is parallel to the unit vector ni, the 
points F, X and M are qr = 0, qx = (^,0) and qM = 
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FIG. 9: For triangular lattice the eigenfrequencies (in units of 
\J D Im) (q) (solid) and CJ2 (q) (dashed) along lines between 
points r, K and M. These points are in the first Brillouin zone 
shown in the inset. 

In Fig. ^1 a square lattice is shown. We now assume 
first and second nearest neighbor interactions (solid and 
dashed lines in the figure) represented by springs of force 
constant D\ and Z?2, respectively. Similar way as before, 



D, 



D, 



Ax 



n, 



^1. 



FIG. 10: Square lattice with first and second nearest interac- 
tions. The straight lines correspond to unstrained springs of 
force constant D\ (solid lines) and D2 (dashed lines) connect- 
ing first and second nearest neighbor atoms in equilibrium 
state. The distance between nearest neighbors is a. The unit 
vectors ni,n2,li and I2 correspond to the directions of the 
unstrained springs. 

one finds that the 2x2 dynamical matrix is 



D(q) = 3- [1 - cos(qnia)] rii (g) 

z— 1 
2 

+ [l-cos(^ql,\/2a^ 



2 V) 

^ - ■ ' ^ l^^l^^ (33) 

where the unit vectors and 1; (i — 1,2) are shown in 
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Fig. El Finally, the two eigenfrequencies can be obtained 
again from Eq. (|^ . In Fig.lTTIthe eigenfrequencies uji (a) 
and UJ2 (q) are plotted along the lines in the first Brillouin 
zone joining the points T, X and M shown in the figure. 
In the Cartesian coordinate system in which is parallel 
to the unit vector ni, the points T, X and M are qr — 0, 
qx = (^,0) and qM = , f ), respectively. 




r X M r 



FIG. 11: For square lattice the eigenfrequencies (in units of 
\fDZJrri] '^i(q) (solid) and ii;2(q) (dashed) along lines be- 
tween points r, X and M. These points are in the first Bril- 
louin zone shown on the right of the figure. The ratio of the 
force constants is D\/D2 = 2. 



VII. CONCLUSIONS 

We first calculated the vibrational modes of the hon- 
eycomb lattice in the harmonic approximation. It was 
assumed that nearest neighbor atoms are connected by 
ideal springs. Using the direct product of vectors we de- 
rived a formula for the spring force acting on an atom. 
The equations of motion for the atoms were then derived 
and the resulting dynamical matrix was given explicitly. 
The vibrational frequencies and modes were determined 
from the eigenvalue problem of the dynamical matrix by 
analytic methods. Our work may provide a starting point 
to the studies of the more complicated lattice dynamics 
of nanotubes. Our general approach is also applied to 
study the lattice dynamics of the equilateral triangular 
and the square lattices. In the latter case, to investigate 



more complicated structures, not only the first but the 
second nearest neighbor interactions are also included. 
Finally, we mention some problems. 

• What is the first Brillouin zone for the square, 
triangular and honeycomb lattice structures dis- 
cussed before? Study the symmetry properties of 
the eigenfrequencies wi(q) and W2(q) in the first 
Brillouin zone. 

• Study the long wavelength limit, that is when 
|q|a ^ 1 for the square, triangular and honeycomb 
lattice structures. What are the eigenmodes in this 
case? Determine the sound velocity. 

• In Fig. El a and b two ladder-type lattices are 
shown. Using the method developed in Sections Hll 
and mil find the vibrational frequencies and modes 
for these lattices. Determine how the results de- 
pend on a in Fig. 1121 b. 

T*^ r*^- Di 

I >< >< I X. """""" °2 

a) 

— ^1 — f — ^ 

/\ :\ :\ :\ Di 

• U ^ * ^ ^ 

b) 

FIG. 12: Square (a) and triangular (b) ladder lattices. Dif- 
ferent types of lines represent different spring constants. In 
the triangular ladder all springs of force constant D2 have the 
same unstrained length. 
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